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ABSTRACT 

We report the results of J-band infrared spectroscopy of a sample of 13 0 = 1 
field galaxies drawn from the Canada-France Redshift Survey, targeting galaxies whose 
redshifts place the rest frame Ha line emission from HII regions in between the bright 
night sky OH lines. As a result we detect emission down to a flux limit of ~ 10“^® ergs 
cm“^ s“^ corresponding to a luminosity limit of ~ 10"^^ ergs at this redshift for a Hq = 
50 km s“^ Mpc,“^ qo = 0.5 cosmology. From these luminosities we derive estimates of 
the star-formation rates in these galaxies which are independent of previous estimates 
based upon their rest-frame ultraviolet (2800A) luminosity. The mean star-formation 
rate at z = 1, from this sample, is found to be at least three times as high as the 
ultraviolet estimates. The dust extinction in these galaxies is inferred to be moderate, 
for standard extinction laws, with a typical Ay = 0.5-1.0 mags, comparable to local 
field galaxies. This suggests that the bulk of star-formation is not heavily obscured, 
unless one uses greyer extinction laws. 

Star-forming galaxies have the bluest colours and a preponderance of dis¬ 
turbed/interacting morphologies. We also investigate the effects of particular star- 
formation histories, in particular the role of bursts vs continuous star-formation in 
changing the detailed distribution of UV to Ha emission. Generally we find that mod¬ 
els dominated by short, overlapping, bursts at typically 0.2 Gyr intervals provide a 
better model for the data than a constant rate of star-formation. The star-formation 
history of the Universe from Balmer lines is compiled and found to be typically 2-3 x 
higher than that inferred from the UV at all redshifts. It can not yet be clearly estab¬ 
lished whether the star-formation rate falls off or remains constant at high-redshift. 

Key words: surveys - cosmology: observations - galaxies: evolution - galaxies: star- 
burst - stars: formation 


1 INTRODUCTION 

The topic of the history of star-formation in the Universe 
has excited much interest in recent years, stimulated by 
the first observations of nearly-normal star-forming galax¬ 
ies at 2 > 3 (Steidel et al. 1996). Previously high-redshift 
studies were limited to highly active galaxies that may be 
poor tracers of the typical star-formation history of the uni¬ 
verse as a whole. Steidel et al. used the colour signature of 
the Lyman break/Lyman-a forest discontinuity being red- 
shifted through optical filters to select high-redshift objects 
(Guhathakurta et al. 1990), these were subsequently con¬ 
firmed spectroscopically on the 10m W.M. Keck telescope. 


Comparison of these objects with the low-redshift {z < 
1) samples of field galaxies (Lilly et al. 1995, Ellis et al., 
1996) appears to show a rise in the Universal star-formation 
rate from 2 = 0 to 2 = 1 and a drop-off at 2 > 3 indicat¬ 
ing a star-formation peak in the 2 = 1-2 epoch (Madau 
et al. 1996). This is also inferred to be the epoch when 
large galaxies with classical elliptical and spiral morpholo¬ 
gies are assembled: Hubble Space Telescope observations in¬ 
dicate they are extant at 2 = 1 (Brinchmann et al. 1998, 
Lilly et al. 1998) but absent in the 2 > 3 sample (Giavalisco 
et al. 1996, Lowenthal et al. 1997). Theoretical developments 
using galaxy formation simulations constrained by the ob- 
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served evolution in the density of neutral gas from Lyman-a 
QSO absorbers show qualitative agreement with this picture 
(Fall et al. 1996). 

However these measurements of star-formation rate are 
based upon the measurement of ultraviolet continuum lumi¬ 
nosity, 1500-2800A in the rest-frame, assumed to be from 
young stellar populations. If dust extinction played a signif¬ 
icant role in obscuring UV radiation they could be under¬ 
estimated by large factors which may change the picture 
completely. 

A more robust way to measure the star-formation rate 
of high-redshift galaxies would be to measure their luminosi¬ 
ties in Balmer recombination lines. This radiation come from 
reprocessed ionising radiation emitted by young stars. This 
approach has two advantages: firstly the ionising radiation 
comes from more massive short lived stars than the softer 
1500-2800A UV and hence falls quickly to zero only 20 Myr 
after star-formation stops. Thus the Balmer luminosity is 
a more direct measure of the instantaneous star-formation 
rate. This contrasts with the UV which continues to rise as 
the stellar populations evolve, typically doubling for exam¬ 
ple between 10 and 1000 Myr at ISOOA. For Ha the main 
dependence is directly on the Initial Mass Function and neg¬ 
ligibly on the temporal evolution. 

The second main advantage is that the Balmer radia¬ 
tion is emitted in the red part of the optical spectrum and 
is thus much less affected by any dust extinction or attenu¬ 
ation than the ultraviolet. For example for typical SMC and 
Milky Way extinction laws (Pei 1992) the ISOOA and 2800A 
extinctions (in magnitudes) range from 2-7 times greater 
than that at Ha (6563A). 

However to observe Ha at high-redshift requires in¬ 
frared spectroscopy, which has not been possible until re¬ 
cently because of the faintness of the sources involved. Pet- 
tini et al. (1998) have secured the first IR spectra of 5 of the 
z > 3 Steidel et al. galaxies, and obtained H/3 luminosities. 
In this paper we report the results of the first measurements 
of the Ha line in a sample of normal 2 ~ 1 field galaxies 
drawn from the Canada-France Redshift Survey (Lilly et al. 
1995, Le Fevre et al. 1995, Lilly et al. 1995B, Hammer et al. 
1995). This sample (hereafter ‘CFRS’) is a highly-complete 
redshift survey of a magnitude-selected {Iab < 22.5) sam¬ 
ple of normal field galaxies. The median redshift is 0.6, and 
galaxies extend out to 2 = 1.3. Because the sample is mag¬ 
nitude selected the 2 > 1 end is dominated by luminous 
L ^ L* galaxies. 


2 OBSERVATIONS AND DATA REDUCTION 

The observations were carried out on May 10-11 and Octo¬ 
ber 3-5 1996 at the UK Infrared Telescope in Hawaii using 
the CGS4 spectrograph (Wright 1994). We chose galaxies 
in the redshift range 0.790-1.048 so that the Ha line would 
lie in the relatively clean part of the near-infrared J-band 
(1.17-1.34 jam) where the atmospheric extinction is relatively 
low (< 15%). We also selected galaxies with detectable [OH] 
emission in the optical CFRS spectra which make up 85% 
of the CFRS sample at 2 ~ 1. 

As well as absorption, the J-band is contaminated by 
numerous airglow OH emission lines, which increase the 
broad-band sky-brightness by a factor of 30 and hinders 


the detection of faint objects. Our observational strategy 
was to observe at high-resolution with CGS4, thus resolv¬ 
ing out the OH background. Because we already knew the 
galaxies’ redshifts from the optical spectra the limited wave¬ 
length coverage at high-resolution was not a problem. More¬ 
over we could exclude galaxies whose redshifts would put the 
Ha line on or close to an OH line. (A similar strategy was 
also adopted by Pettini et al.) We observed with the 150 
lines/mm grating and the 1 arcsec/1 pixel slit giving a res¬ 
olution A/AA = 2200. We determined empirically that at 
this resolution OH lines contaminated 50% of the bandpass, 
i. e. we had to exclude 50% of the high-redshift galaxies, but 
in the remaining clean part of the bandpass the mean back¬ 
ground was only 20% that of broad J. 

Targets were acquired using the following procedure: 
first we ‘peaked-up’ on a very bright star (U = 1 — 2 mags) 
within 1-2 degrees of the target. This involves centering the 
star on the optical finder TV, reading the IR array in a 
continuous ‘MOVIE’ mode, and then adjusting the offset 
between the axis of CGS4 and the TV until the IR flux 
is maximised. This assures the IR slit is aligned with the 
TV crosshair. Then we went to a fainter star, typically 17th 
mag, within 1-2 arcmin of our target galaxy and measured 
off the same coordinate system, centered the star on the 
TV crosshair and did a blind offset on to the target galaxy. 
(The targets were too faint to see on the TV). We would 
then autoguide either on the offset star or another bright 
star in the region. 

Observations were made stepping the InSb detector ar¬ 
ray in 0.5 pixel increments to fully sample the instrument 
profile and nodding the telescope ± 9 arcsec along the slit 
between ‘OBJECT’ and ‘SKY’ positions to facilitate sky- 
subtraction (though note the object is still on the slit in both 
positions). Individual exposures ranged from 10-15 minutes. 
The typical seeing was 1.0 arcsec. A total of 13 objects were 
observed: these are listed with their total exposure times in 
Table Standard wavelength calibration and Hatfield cor¬ 
rections were applied. The October observing run was af¬ 
fected by the spectrograph slit being jammed out of position 
which caused the lines to be tilted on the image. The shifts 
were measured by cross-correlation and the tilt corrected by 
r e- int erp olat ion. 

Even with the resolved OH background accurate sky- 
subtraction is critical to detection of faint lines. To first or¬ 
der the sky can be removed by simply subtracted the pairs 
of offset frames, though this leaves residual signal due to 
temporal sky changes. To second order the residual sky was 
removed by performing a polynomial interpolation along the 
slit, excluding the two object regions, and subtracting. This 
leaves no systematic residual, though the regions near the 
OH lines are still noisier due to the extra Poisson contri¬ 
bution, with the result being a 2D image with a positive 
spectrum in the ‘OBJECT’ row and a negative spectrum in 
the ‘SKY’ row. For each image we also made a pixel mask 
to exclude bad pixels on the detector and regions with noisy 
OH residuals. 

In many of the images there were strong Ha and con¬ 
tinuum detections. We summed all the Ha lines with good 
detections, fitted Gaussians spectrally and spatially to de¬ 
fine the typical line profile and then used this mean profile 
(with the pixel mask) to optimally extract all bright, faint 
and possibly non-detected objects in a consistent manner. 


© 1998 RAS, MNRAS 000, Mi 


Ha star-formation rate at z = 1 3 



1.335 1.34 

X / microns 



1.19 1.195 1.2 

X / microns 




1.19 1.195 1.2 


X / microns 


'■ 23.1313 1- 

M DET 

1 ^ : 



1.19 1.195 1.2 


X / microns 



1.24 1.245 1.25 


X / microns 



1.25 1.255 1.26 


X / microns 




1.34 1.345 1.35 


X / microns 




1.34 1.345 1.35 


X / microns 



X / microns 



1.185 1.19 1.195 


X / microns 

Figure 1. The J-band spectra of our CFRS galaxies. The spectra have been continuum subtracted and centered on the Ho line. The 
large central arrow indicates the predicted Ho position based on the optical redshift, the associated horizontal error bar denoting the 
1(7 error on the redshift. The two smaller arrows show the predicted positions of the [Nil] 6548A and 6583A lines. The vertical dotted 
lines indicate the positions of strong night sky OH lines which have been masked out. The objects where we find a > 2(7 Ho detection 
are labelled ‘DET.’ 


In most of our observations we found that the nega¬ 
tive spectrum was typically weaker, or even absent, com¬ 
pared to the positive spectrum. This is attributed to the 
fact that any small errors in acquisition are magnified when 
stepping 9 arcsec away from the centre along the slit as the 
rotation is not precisely known. When there was a signifi¬ 
cant negative spectrum it was combined in a weighted man¬ 
ner (to maintain consistent exposure times) into the pos¬ 
itive spectrum. This marginally increases the signal/noise, 


though none of the flux measurements presented below are 
significantly changed if this step is omitted. 

Next the spectra had atmospheric absorption removed 
using a smooth-spectrum standard and were flux calibrated 
using flux standards. 

Finally we applied aperture corrections to allow for our 
finite 1.0 arcsec wide slit. Our objects are resolved, our mean 
Gaussian spatial profile along the slit of the Ha line has a 
FWHM of 2.0 arcsec (which is consistent with the typical 3- 
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5 arcsec isophotal optical diameters measured for the CFRS 
sample (Hammer et al. 1997)). The flux calibration stars are 
observed through the same slit in 1.0 arcsec seeing. Assum¬ 
ing Gaussian profiles for both we derive a relative correction 
factor of 1.7, by which we multiply our spectra. While this 
can only be a rough correction, it gives fluxes consistent with 
the optical lines (see section ^). It is also only of order unity, 
so even if ignored our conclusions below are not drastically 
altered. 

The continuum level was determined in each case by 
taking the mean of the data points ± 0.02 microns either 
side of the emission line, ignoring the masked points and 
the emission line. The noise level of the data was determined 
empirically from the RMS in the same region. In most cases 
we found significant continuum emission. Note that for indi¬ 
vidual pixels the continuum is mostly below the noise, it is 
only by summing up we get a significant detection. To check 
the validity of our measurements we repeated the same pro¬ 
cedure for an off-object row in the longslit spectra. This 
give non-detections in all cases, so we are confident in our 
procedure. 

After subtraction of the continuum level from the data 
we computed the line flux by summing the flux ±N pixels 
around the line, excluding the masked pixels, where N = 7 is 
chosen to be the typical line FWHM. This is done regardless 
of whether there appears to be a detection or not (as the 
expected wavelength is known). 

The final line fluxes and continuum fluxes (converted 
to AB mags and denoted Jc) are given in Table along 
with the most useful CFRS parameters of our objects. Ha 
fluxes < Icr are set = 0. We find 7 detections > 2a and 5 
detections > 3a out of our 13 objects. 

We convert these to luminosities using a Hq = 50 km 
s“^ Mpc,“^ and go = 0.5 cosmology in Table We use this 
cosmology for the rest of our paper. We note that in our 
further analysis our conclusions are based on comparing the 
luminosities of the same galaxies at different wavelengths, 
since all the galaxies in our sample lie close to 2 = 1 our 
conclusions are essentially unchanged if we use a different 
cosmology. 

To complete the quantities derived from the lines we fit 
Gaussians profiles (excluding masked pixels) to derive ve¬ 
locity line widths for all our detections. Each of the fits was 
checked visually by plotting on top of the line; the instru¬ 
mental resolution was determined by fitting to unresolved 
night-sky OH lines in the region of the galaxy line and the 
value was subtracted in quadrature from the galaxy line 
widths. It should be noted that all our lines, but one, are 
well-resolved as we expected given the spectral resolution 
and the typical velocity widths of galaxies. The spectra line 
Full Width Half Maximum (FWHM) of the galaxies range 
from 3-5 pixels (one pixel = 2-3 A depending on the spec¬ 
trum) compared to 2.4-2.6 pixels for the sky lines. The re¬ 
sults are presented in Table we make no detailed analysis 
here, we just note that the typical velocity FHWMs are in 
the range 200-400 km/sec expected for large L ~ L, galax¬ 
ies. 


O 



Figure 2. Comparison of line luminosities of the sample in [Oil] 
and Ho. Dotted lines show slopes of Ho/[OII] = 1, 2, 3, the spread 
of ratios found be Kennicutt (1992). 


3 COMPLETENESS 

In 8/13 of the galaxy spectra Ha emission was detected; 
additionally 9/13 of the galaxies had detected continuum 
emission. This proves quite useful: we can check if our line 
non-detections are due to poor acquisition by comparing the 
continuum level of our line detections and non-detections. 
There are only two cases in which there is no line and no 
continuum detection. The general trend is that the average 
continuum flux is brighter for the non-detections, this trans¬ 
lates in to an median {I—Jc)ab colour 1.1 mags redder. This 
argues against the slit missing the object, in fact the trend 
towards redder colours is precisely what is expected for non¬ 
line emitting objects. We note that at 2 = 1 a Scd galaxy 
should have observed colours (7 — J)ab ~ 1.0 and an E/SO 
galaxy (7 — J)ab = 2.0, using template SEDs from Ken¬ 
nicutt (1992). We find median colours of (7 — J)ab = 1.6 
for the Ha detections and (7 — J)ab = 2.7 for the non¬ 
detections, which agree very well given our actual galaxies 
will differ in detail from Kennicutt’s templates. 

Figure ^ shows the final set of spectra, which have been 
continuum subtracted and centered on the Ha line. The line 
is in all cases found within the error box given by the optical 
redshift (A 2 ~ 0.002). Our resolution is high enough that 
Ha is well separated from the [Nil] lines so we do not have to 
correct for blending. Also in a few cases (e.g objects 22.1406 
and 22.070) we see evidence for one of the weaker [Nil] lines 
as well as the main Ha line. This is additional evidence for 
the robustness of our detections. Note in many cases the 
position of one or both of the [Nil] lines is occluded by a 
noisy OH residual. 

We can also test our completeness using GFRS val¬ 
ues for the [OH] flux (Hammer et al. 1997) to calculate a 
Ha/[OH] line ratio. In the GFRS sample in this redshift 
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range 85% of galaxies have [Oil] emission. Hammer et al. 
tabulate the [Oil] equivalent width and flux, the latter is 
aperture-corrected by comparing the spectra at 5500A 
with their H-band image photometry. Other than exclud¬ 
ing galaxies with zero [Oil] emission we made no attempt 
to concentrate on objects with the strongest [Oil] emission. 
Thus the mean [Oil] equivalent width (32A) and range (10- 
60A) of our small sub-sample are consistent with random 
sampling from the larger sample. The CFRS galaxies at 
z = 1 with zero [Oil] are all at the red end of the colour 
distributions mV — I and I — K. Thus the colours also in¬ 
dicate they are not signihcant star-forming systems and we 
conclude their exclusion has no significant impact on our 
conclusions below. 

One might also ask the question: is the lack of Ha de¬ 
tection in some of our objects consistent with the presence 
of [OH] in our optical spectra? This is addressed in Figure 
^ where we plot the strength of the two lines against each 
other. It can be seen that given the Ha error bars all points 
are consistent with a reasonable linear correlation. Kenni- 
cutt (1992) estimates the line ratio Hq/[OH] as having a 
median value of about 2 (with a spread from 1 to 3) for 
a sample of local star-forming galaxies (1.0 mag mean ex¬ 
tinction). Our observed values at high-redshift are entirely 
consistent with Kennicutt’s median and spread, implying we 
too have small amounts of extinction which agrees with our 
hndings below in Section The points with zero Ha are 
consistent with detected [OH] given the larger error bars. 
Note the line ratios are also good evidence that our aper¬ 
ture corrections are reasonable, if omitted we would obtain 
a much lower ratio Ha/[OH] = 1. 


4 STAR FORMATION RATES FROM Ha AND 
UV 

Using models of population synthesis it is possible to cal¬ 
culate the relation between input star-formation rates and 
output UV and Ha fluxes. The basic principle is that the 
UV light is dominated by short-lived main sequence stars 
(the Ha light is reprocessed ionising UV) so the number of 
them in a galaxy is proportional to the star-formation rate. 
The prescription for this calculation is simple. For reference 
we outline it (and the corresponding assumptions) in detail: 

(i) For a given time-dependent star-formation rate the 
population synthesis code gives the UV stellar spectrum as 
a function of time. Note that Kennicutt (1983) in deriving 
his calibration uses a simple grid of stars of different masses 
with evolutionary tracks current at the time. 

(ii) For the UV continuum estimators one takes an aver¬ 
aged flux (e.g. through a synthetic box filter) at a specific 
wavelength (e.g. 1500A, 2800A). At the longer wavelengths, 
with increasing stellar lifetimes, the conversion from a con¬ 
stant star-formation rate to a UV flux is age-dependent (see 
for example Pettini et al. who uses different factors at 1500A 
for 10^ and 10® years). One then also needs to apply correc¬ 
tions for dust attenuation and/or Lyman absorption by line 
of sight systems (e.g. Madau 1995). 

(iii) For Ha one calculates the number of ionising Lyman 
continuum photons. This flux comes from the most short 
lived stars radiating at A < 912A with lifetimes of < 10 
Myr. Then it is assumed all this radiation is absorbed by 
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Figure 3. Illustration of the dependence of the UV spectrum 
on metallicity and time for a constant star-formation rate. These 
are theoretical galaxy spectra from the BC96 population synthe¬ 
sis models (Salpeter IMF). The position of the 1500A and 2800A 
±10% bands are marked along with the ionising photons which 
are reprocessed into Ha. The vertical order of the labels is the 
same as the vertical order of the spectra. Panel (a) shows the 
spectrum at 1 Gyr for a decreasing metallicity sequence. Panel 
(b) shows the solar metallicity spectrum for an increasing time 
sequence. It is easily seen that while the redder UV bands are 
relatively insensitive to metallicity they are strongly time depen¬ 
dent until 1 Gyr old. The reverse is true for the ionising photons 
— the stellar spectra below 912A are identical for times > 20 
Myr at a given metallicity. Hence the resulting Ha flux is mainly 
dependent on metallicty. 


intervening hydrogen gas in the galaxy in which the forming 
stars are embedded and that none leaks out. In practice it 
seems this is very close to the truth. At very high-redshift 
the Lyman limit can be observed in the optical (e.g. the 
Steidel et al. galaxies) and the flux does indeed go to zero, 
but as these galaxies are identified by the Lyman breaks 
this could be a selection effect. For local galaxies there have 
been limited A < 912A spectroscopy, for example Leitherer 
et al. (1995) observed a sample of 4 starburst galaxies with 
the Hopkins Ultraviolet Telescope and concluded that < 3% 
of the ionising photons escaped. Constrained models of the 
ionising radiation field of the Milky Way (Bland-Hawthorn 
& Maloney 1997, 1998) indicate that approximately 5% may 
escape. 

(iv) The ionising photons are reprocessed into recombi- 
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nation lines, and the relative strengths can be calculated 
in detail. Hummer & Storey (1987) calculate 0.45 Ha pho¬ 
tons are emitted per Lyman continuum photon for case B 
recombination. This number is quite robust, over a range of 
nebulosity conditions (10^ —10®K, 10^ — 10“^ electrons cm“^). 

As there are a number of values for these conversion fac¬ 
tors in the literature we thought it would be useful for ref¬ 
erence purposes to systematically tabulate these for a set of 
models. This also serves to illustrate the range of variations 
and trends. The results of our calculations are given in Ta¬ 
ble ^for UV 1500A, 2800A and Ha conversions. The Bruzual 
& Chariot (1993, 1996) models (‘BC96’) offer a range of 
metallicities (albeit using theoretical model atmospheres — 
the ‘kl96’ models). The PEGASE models (Fioc et al. 1997) 
offer two sets of post-main sequence evolutionary tracks but 
only solar metallicity. Both sets of models offer several Stel¬ 
lar Initial Mass Functions (IMFs); we tabulate the results 
using the IMFs of Salpeter (1955) and Scalo (1986) because 
these bracket the upper and lower limits of the variation. 
We also tabulate some other values given previously in the 
literature. 

For all these models Ha reaches a constant asymptote 
after 20 Myr, then stays constant to within a 5-10% out 
to 10 Gyr. Thus it can be used as a tracer of the instanta¬ 
neous star-formation rate on time-scales of order ~10 Myr. 
For reference we give the value at 100 Myr. The conver¬ 
sion is very sensitive to the IMF, as the ionising flux comes 
from the most massive stars, with Salpeter giving a value 
2.8-3.4 X higher than Scalo for the different models. An¬ 
other important point, which has not been remarked upon 
in the literature discussing high-z galaxies, is the strong ef¬ 
fect of metallicity below 912A. This occurs in the stellar 
evolutionary models because the most luminous stars (con¬ 
tributing below 912 A), have a much greater dependence of 
effective temperature on metallicity than the less luminous 
stars contributing at longer wavelengths. The effect of the 
hotter stars at low metallicty is to raise the Ha flux about 
1.7 X compared to the solar metallicity value for a given 
star-formation rate. 

The average UV luminosity per unit Hz is computed as 
the total energy integrated through a ± 10% box filter di¬ 
vided by the total bandwidth. This gives a negligible colour 
term. For a 10% filter the difference in the average between 
a constant fi, spectrum (which approximates these spectra 
for > 912 A) and a much redder constant f\ spectrum is 
only 1%. The IMF dependence is somewhat less than for 
Ha (typically a factor of 1.5-2.5) as these UV bands come 
from less massive stars and the metallicity dependence on 
temperature is much smaller. 

However there is a stronger time dependence, with the 
flux increasing significantly between 10^ and 10^° years, es¬ 
pecially at 2800A and for the Scalo IMF which is much less 
rich in massive stars. For cosmological times (> 1 Gyr) 
the I 500 A flux is relatively insensitive to time as noted 
by Madau et al. (1998) (except for the extreme 0.02 So¬ 
lar metallicity case). Thus it can be used as a tracer of the 
star-formation rate on time-scales of 1 Gyr. The 2800A flux 
is not as stable and changes at the 10-40% level over 1-10 
Gyr. There are further < 50% variations in the exact con¬ 
version depending on the model and metallicity assumed. 


o 



L(2800A) / 10^’^ergs Hz“^ 

Figure 4. Comparison of the Ho vs UV continuum flux at 2800 A 
for the individual galaxies. The overlayed axes show the locus of 
constant star-formation rate for a set of fiducial models cover¬ 
ing the range of UV/Ho ratios given in the last column of Ta¬ 
ble ^ The numeric labels on the axes are the corresponding star- 
formation rates in MQyr~^. The dotted line shows the ratio of 
the mean luminosity densities derived in Section 

These time and metallicity effects are illustrated graph¬ 
ically in Figure Finally it is worth noting that the range 
in the ratio of Ha to UV is less than the absolute range; as 
both come from high-mass stars the choice of IMF matters 
somewhat less. 


5 THE STAR-FORMATION RATE AT Z = 1 

Turning back to our data we can do a direct galaxy by galaxy 
comparison of the star-formation rates inferred from UV and 
Ha. 

Our UV fluxes are particularly robust because for red- 
shifts near unity the rest frame 2800A light corresponds very 
closely to the observed frame U-band light. To a first approx¬ 
imation we can simply ignore the K-correction and derive 
the 2800A flux directly from the CFRS U-magnitudes. We 
refine this slightly by using the SED fits from Lilly et al. 
(1996); this corrects the fluxes by 10-20%. Thus the UV lu¬ 
minosities we use are the same as the raw data which goes 
into the 2 = 1 star-formation rate determinations of Madau 
et al. (for 0.2 < 2 < 1 this was based on the GFRS luminos¬ 
ity density functions of Lilly et al. 1996) We use a formal 
error bar of 20% for our UV fluxes to give errors represen¬ 
tative of the V photometric errors and the K-corrections. 

We plot the Ha vs UV luminosities in Figure ^ and 
overlay lines for a set of conversions from Table Note these 
are only valid for a continuous, constant star-formation rate 
— the more complex problem of bursts is considered below 
in Section ^ 

It is clear from the plot that there is a order-of- 
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magnitude agreement in the star-formation rates derived 
from the two methods and a reasonable correlation, i.e. 
strong UV systems are usually strong Ha systems. There 
are no galaxies with an extremely large excess of Ha rel¬ 
ative to UV, which would occur if the star-formation was 
highly obscured by dust. 

We do Hnd points with zero Ha flux and appreciable UV 
flux, which we would expect to detect within the measure¬ 
ment errors if there was a linear correlation. This can not 
be due to dust as the latter can only enhance the UV/Ha 
ratio not diminish it. As noted in Section ^ these spectra 
are detected, they just do not contain significant Ha. One 
physical effect that will explain this is the relative lifetimes 
of stars contributing to the UV and Ha as mentioned in Sec¬ 
tion When one moves away from the simple scenario of 
constant star-formation rate the picture changes consider¬ 
ably. For example for a instantaneous starburst the Ha flux 
will drop to effectively zero after 30 Myr while signihcant 
UV flux will persist up to 1000 Myr. This will produce low 
Ha points such as are seen in our sample. Weak Ha should 
be present in all our galaxies because they are known to have 
[OH] emission (see Section ^, however the UV/Ha ratio will 
scatter more widely and it is certainly possible for Ha to be 
below our detection limit despite significant [OH] and UV 
emission. To quantify these effects requires us to develop a 
proper model for starburst activity in galaxies. We do this 
below in Section 

Nevertheless when one averages over many galaxies the 
effect of bursts should cancel out in the mean — the con¬ 
tinuous star-formation conversions of Table ^ should be ap¬ 
plicable for an ensemble of galaxies. (This approximation is 
tested below using the methods developed in Section ^ and 
is found to accurate to ~ ±10%, even for this small sam¬ 
ple.) It can be seen that the values scatter about a ratio of 
2-3 times as much star-formation inferred from the Ha as 
UV. This is balanced to some extent by the zero points. By 
summing over all the galaxies we find the following relation 
between the luminosity per comoving volume in Ha and UV: 


L(2800A) 

L(Ha) 


3.1 ±0.4 


in the units as in the last column of Table ^ This ratio is 
plotted as the dotted line in Figure |^. 

The ratio is somewhat lower that those predicted by 
the models which give values in the range 4-14 for 0.2-2.5 
Solar metallicities. What can cause this discrepancy? It is 
too large to be encompassed by our range of metallicities; 
one can boost the UV/Ha ratio if we assume a younger 
age than ~ 1 Gyr. However changing to 0.1 Gyr only in¬ 
creases the slopes of the solid lines in Figure ^ by ~ 30%; 
we could achieve the observed slope if star-formation is only 
^ 0.01 Gyr old in all the detected galaxies but this would 
be unlikely simply due to random sampling — the redshift 
range of the sample corresponds to a timespan of ~ 1 Gyr 
at 2 = 1. The most likely explanation for the discrepancy is 
the presence of dust attenuating the ultraviolet light. 

If we adopt fiducial model values of (10.9, 6.3) for 
(Scalo, Salpeter) for Solar metallicity, which agree well be¬ 
tween PEGASE and BC96 models (kl96 tracks), the ratio 
of Ha to UV inferred star-formation rates are (3.5, 2.0). Us¬ 
ing Pei (1992) extinction formulae (a standard dust-screen 


model) we then derive a mean sample Av of (1.0, 0.5) mags 
for the SMC law, (1.1,0.6) for the Galactic law. 

These extinction values are entirely consistent with 
what is found from studies of star-formation from Balmer 
lines in local normal Sab galaxies. For example Kennicutt 
(1983) found Ay = 1.0 mags and a study of low redshift 
2 < 0.3 CFRS galaxies using the ///3/Ha line ratios by 
Tresse and Maddox (1998) found the same value. 

Calzetti et al. (1994) and Calzetti (1997) have proposed 
an empirical dust-attenuation law for heavily reddened star- 
bursts {Av — 2.2 mags). In this model the nebular lines 
have about twice the optical depth as the stellar contin¬ 
uum; moreover although the lines are well described by a 
standard Galactic screen law the stellar continuum is em¬ 
pirically described by a greyer law in which the dust and 
stars are intermixed. Using the Calzetti law we derive Ay 
attenuations of (2.6, 1.4) mags (for stars; xl.4 for nebulae); 
because the attenutations of the Ha and 2800A continuum 
are more similar than for the screen models a much higher 
obscuration is required to match the observed excesses; at¬ 
tenuations of (2.9, 1.6) mags result for Ha and (4.2, 2.3) 
mags for 2800A stellar continuum. It should be noted how¬ 
ever that it is still not established that corrections derived 
for the Balmer lines in starburst regions of nearby galaxies 
are appropriate to the integrated light of the distant galaxies 
studied here. This question remains open. 

We can now examine the total correction for dust in 
both the Ha and UV determined star-formation rates. The 
Milky way law gives corrections for (Scalo, Salpeter) IMFs of 
x(2.2,1.5) for Ha and (8.0,3.1) for the UV. These numbers 
are the same to within ±10% for the SMC extinction law; 
this is because the two extinction curves differ most strongly 
in the UV at << 2800A (e.g. at the 2175A dust feature) 
and in the near-UV and optical they are very similar. The 
Calzetti law of course gives much larger values: the final 
star-formation rates are an additional factor of ~ 6x higher 
for the Scalo IMF and ~ 3x higher for the Salpeter IMF. 

Finally we can compare the star-formation rate of our 
2^1 CFRS galaxies with local counterparts. We adopt the 
Salpeter IMF as that is conventionally used for deriving the 
local rates. For the range of dust corrections we have derived 
we obtain rates of ~ 20-60 Mq yr“^, comparable to local 
starbursts (e.g. Calzetti 1997) and much greater than the 
typical 4 Mq yr”*^ found for local normal spirals (Kennicutt 
1983) and the Milky Way (Smith et al. 1978). 


6 MODELING OF STARBURSTS 

As noted earlier the interpretation of Ha and UV luminosi¬ 
ties as star-formation rates becomes more complicated if 
non-constant star-formation histories are assumed. For this 
reason we developed a mathematical frame-work for explor¬ 
ing this and to see how well we could reproduce the observed 
distribution. 


6.1 Methods 

The principles are based upon maximum likelihood and are 
a 2D generalisation of the methods developed by Abra¬ 
ham et al. 1999 for colour-colour fitting. For a given star- 
formation history we can run a spectral evolution code and 
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Figure 5. Results of fitting burst models to our extinction-corrected Ho and UV (2800A) data. The three rows show three different 
models as described in the text, they are: continuous star-formation (top row), fixed mass exponential bursts plus continuous star- 
formation (middle row) and variable mass exponential bursts plus continuous star-formation (bottom row). The left hand column of 
panels shows the observational data (points) compared with a model distribution generated from the best fit model (contours correspond 
to a factor of 10 in probability density). The middle column of panels show the likelihood contours of the main burst parameters generated 
by our fitting. (Contours correspond to a factor of 10 in likelihood, the circle marks the maximum likelihood point.) The right hand 
column show a realisation of the best fit model, i.e. 13 simulated data points generated from the model distribution with the observed 
errors (points < Icr are set = 0 as in the data). Note the continuous star-formation includes a component of galaxy—galaxy scatter. 


calculate Ha and UV luminosity as a function of time using 
the methods of Section We evolve the models for 5 Gyr 
and sample at random times (we can do this because all the 
galaxies are at similar redshift) to create a 2D model dis¬ 
tribution of light in the the (Ha, UV) plane. The code then 
computes the likelihood of drawing the observational data 
from the model: 




r,,u N ( {hi-hf 
P{h,n)e.p I 


' —oo —oo 


2nAhiAui 


{Ui - 

2 Auf J d.u dh 


where the observational (Ha, UV) points are {hi,Ui) 
with errors {Ahi, Aui), and the model points have proba¬ 
bility density P{h,u). 

The star-formation histories are parameterised and like¬ 
lihood space is explored to find the maximum likelihood. We 
use an adaptive algorithm where a large parameter space is 


explored with a coarse parameter grid to locate the peak 
and the region around the peak is then examined with a 
finer grid to give confidence limits on the fitted parameters 
from AC. 

Once we have the best fit parameters for a model, we 
can then create simulated observational datasets. This is 
done 1000 times, the maximum likelihood fit being recom¬ 
puted each time, this allows us to normalise our relative 
likelihood into an absolute probability of the observed data 
given the model. 


6.2 Star-formation histories explored 

We parameterised the star-formation histories as continuous 
star-formation, with galaxy to galaxy scatter, plus a ran¬ 
dom distribution of bursts. Initially the total star-formation 
rate is kept constant and normalised to the values derived 
in Section p. This is a good approximation as an ensemble 
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Figure 6. The star-formation rate, luminosity history in H« and 
UV ( 2800 A) typical of our best fit models. This is for the case 
of exponential bursts of variable mass, with a continuous star- 
formationcomponent, i.e. the model shown in the lower 3 panels 
in Figure ra. It can be seen that for long inter-burst periods galax¬ 
ies are indeed quiescent in H«, but not the UV. 


of galaxies undergoing bursts will approximate continuous 
star-formation and results in one less free parameter. For 
our further analysis we confine ourselves to the BC96 mod¬ 
els (Solar metallicity, ‘kl96’ atmospheres) and the Salpeter 
IMF (as the latter gives better results for fitting the star- 
formation histories, colours and mass/light ratios of galaxies, 
see for example Kennicutt (1983) Madau et al. 1998, Calzetti 
1997, Lilly et al. 1996). We correct all (Ha, UV) points using 
the Av = 0.6 mags Milky Way law derived in Section]^; the 
effect of adopting the Calzetti law is just a simple scaling to 
globally 3x higher star-formation rates and does not affect 
the details of the analysis. 

To start with we modeled a simple continuous star- 
formation model to check our code was giving sensible re¬ 
sults. To make it more realistic we introduced a scatter, Ca, 
between galaxies following a normal distribution (slightly 
corrected to avoid negative star-formation rates). 

The best fit results are shown in Figure ^ as a likelihood 
contour plot overlayed on the model points. For comparison 
we also show a simulated set of data points drawn from 
the model distribution. It can be seen that some scatter is 
introduced; this originates physically from the time variation 
of UV light even for constant star-formation. The best fit 


parameter values are shown in Table ^ — the model is a 
very bad fit to the distribution of data. 

One might ask if the scatter could be explained by vari¬ 
ation in extinction between galaxies. While it is possible 
this can explain some of the variation it can not explain 
the points with large amounts of UV emission and small 
amounts of Ha emission. This is because dust quenches the 
UV much more than the Ha — precisely the opposite effect 
to that sought. We next explored the effects of starbursts to 
see if these could plausibly explain the observed distribution. 

Initially we tried fixed-mass bursts, then we tried a 
scheme for allowing the burst masses to vary. Since we could 
find little information in the literature as to an appropriate 
mass function to adopt for bursts we invented our own sim¬ 
ple phenomenological scheme: bursts are parameterised by 
a mean mass (M^) and standard deviation {Me). The dis¬ 
tribution is assumed to be normal. While this has no phys¬ 
ical basis it only has 2 parameters and at least allows us 
to investigate the effects of mass distributions. and our 
global star-formation rate normalisation fixes the mean in¬ 
terval between bursts. Finally we have F, the fraction of 
the star-formation occurring in bursts. For the form of the 
bursts we consider two cases (modeled after BC96): a con¬ 
stant burst of length r and exponential bursts of e-folding 
time T. 

The results of this exercise is shown in Table ^ It can be 
seen that the burst models provide much better fits than the 
the continuous models. The Monte-Carlo realisations show 
that the latter generate data sets like the observed one only 
about 1% of the time. This is because the continuous star- 
formation can not generate enough variation in the UV/Ha 
ratio. The burst models are much better and generate syn¬ 
thetic points which look like the data ^ 10-20% of the time. 
This is because the variation in the star-formation rate is 
the principle cause of variation in the UV/Ha ratio. Allow¬ 
ing the burst mass to vary improves the fit only slightly, we 
conclude that our data does not constrain the shape of the 
burst mass function significantly. The best fit burst fraction 
is ~ 1, indicating the burst mode is preferred. The results 
are illustrated graphically in the lower two rows of Figure ^ 
which compares the dust-corrected observational data with 
model realisations for a sample of key models. The likelihood 
contours of the main parameters are also shown. 

The best fit mass of a typical burst is 2-5x10® Mq, 
corresponding to a time interval between bursts of typically 
~ 200-300 Myr and the characteristic time r is of ~ 100-200 
Myr. This mean the bursts usually overlap in time. These 
values are similar to what one might expect intuitively based 
on the data: the chance of catching a galaxy in the Ha qui¬ 
escent stage has to be of order 1/3 to reproduce the fraction 
of points seen with UV but no Ha. 

We note that this kind of ‘continuous but episodic’ star- 
formation with several bursts per Gyr is of the same form 
as that found in local starburst galaxies by Calzetti (1997). 
Our average star-formation rates are of the same order too. 
In between bursts the star-formation rate and Ha flux does 
indeed drop close to zero while the UV persists (see Fig¬ 
ure ^ which shows a sample time-dependence) due to the 
stellar lifetime effects mentioned in Section Thus the zero 
Ha points in our data (given the error bars) are naturally 
explained. 

Finally with these tools we were able to test how well 
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Figure 7. Postage stamps (10.1 X 10.1 arcsec) of the six galaxies having images from our HST morphology programme. 


an ensemble of galaxies converged to approximating a con¬ 
tinuous star-formation rate, as assumed in Section H. To do 
this we re-ran the likelihood fitting, this time fitting for the 
total star-formation rate as a free parameter. The results of 
this gave rates of between 19 and 23 Mq yr“^ per galaxy, 
which agrees well with the value of 20 Mq yr”*^ calculated 
in Section ^ for the same Salpeter IMF. 

While these simple models could do with some elabora¬ 
tion to obtain a better fit, we are near the limit of what can 
be inferred from 13 data points. It is clear this sort of de¬ 
tailed approach will benefit greatly from future observations 
and much larger samples. 


7 MORPHOLOGICAL TRENDS 

Six of the galaxies in our sample have morphological in¬ 
formation from our programme of Hubble Space Telescope 
high-resolution imaging of CFRS and LDSS2 high-redshift 
galaxies (Brinchmann et aZ.1998, Lilly et al., 1998). This 
is obviously an even more limited sample, but we can look 
qualitatively at the dependence of star-formation rate on 
galaxy type. 

The classifications are listed in Table 1 and postage 
stamps of the galaxies are shown in Figure 7. There are 3 


galaxies classified as ‘Peculiar’. All have detected Ha emis¬ 
sion and blue colours {{V — I)ab < 1) and one (14.0600) 
has the highest star-formation rate in our sample. Two of 
these are classed as ‘mergers’ and one as a close pair indicat¬ 
ing an association between star-formation and interaction. 
There are two galaxies classed as ’Compact’. Both of these 
are also blue {{V — I)ab < 1.2) and it is interesting to note 
that while one has quiescent Ha and strong UV the other 
has a large Ha excess (3.2). Finally we have 03.0316 a red 
spiral {{V — I)ab = 2.9) which is quiescent in both UV and 
Ha. 

Finally we note that galaxies with star-formation, 
whether inferred through Ha or UV, have the bluest V — I 
colours {{V — T)ab < 1.2). This is not in contradiction with 
the modest extinction values derived in Section — an ex¬ 
tinction of Av = 0.6 mags would only redden the observed 
V — I (rest 2800A — B) by a small 0.4 mags whereas the 
difference between and old and young stellar population is 
of order 2-3 mags. 


8 COMPARISON WITH OTHER RESULTS 

There are now enough measurements of Balmer line star- 
formation rates at high and low redshift to construct the 
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Figure 8. Star-formation history of the Universe inferred from Balmer lines (Ha and H/3) compared to UV and far-IR/sub-mm determi¬ 
nations. The UV points are: Treyer et al. (1998; open triangle), Lilly et al. (1996; open circles), Connolly et al. (1997; open stars), Madau 
et al. (1996, 1998; open squares). The Balmer points are Gallego et al. (1995; filled triangle), Tresse & Maddox (1998; filled diamond), 
this work (filled circle), Pettini et al. (1998) correction of the UV Madau et al. point at z = 2.8 (filled square). We also plot the sub-mm 
derived point of Hughes et al. 1998 (solid star) and the far-IR ISO derived points of Flores et al. (1998) (open cross). 


first star-formation history of the Universe in Balmer light 
to compare with the previous UV measurements. For con¬ 
sistency we use the Salpeter IMF, Ho = 50 kms“^ Mpc“^, 
go = 0.5 throughout. All points are re-derived from their 
original luminosity densities in a consistent manner using 
the UV, Ha factors in Table ^ for BC96 (kl96) with Solar 
metallicity. 

This is shown in Figure The point at 2 = 0 comes 
from the Gallego et al. (1995) local objective prism survey 
and is based on Ha. Tresse and Maddox (1998) have mea¬ 
sured the Ha luminosity function at 2 = 0.2 from the CFRS, 
at which point Ha is still available in the optical CFRS spec¬ 
tra. They find a value a factor of two higher than the UV 
measurements. 

Our Ha measurements are used to derive a new value 
for the star-formation rate density in the CFRS at 2 = 1. 
This is higher by a factor of 3.lx than the UV point. At 
2 = 2.8 we show the point derived from the work of Pettini 


et al., who used CGS4 to measure the H/3 line in 5 of Steidel 
et al.’s galaxies. They infer star-formation rates 0.7-7 times 
higher than derived from the UV at 1500A rest and typical 
extinctions A(1500A) = 1-2 mags. We show this as a factor 
of 3 above the UV point. As well as the Balmer lines we also 
plot the point of Hughes et al. based on sub-mm observa- 
tious and the points of Flores et al. (1998) from far-infrared 
ISO observations. It should be noted though that the deriva¬ 
tion of the latter are qualitatively different from the Balmer 
line and UV measurements: the far-IR and sub-mm bands 
measure UV reprocessed by dust into thermal radiation and 
hence they are sensitive to galaxies which might not appear 
at all in the optical. Moreover the Hughes et al. points are 
based on an assumed redshift distributions for sources which 
have not yet been verified and has been disputed (Richards 
1999). Finally it is also worth noting that the Hughes et al., 
Connolly et al. and Madau et al. points are all based on the 
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same patch of sky, the Hubble Deep Field North (Williams 
et al. 1996), which may not be representative. 

It can be however seen that the general trend is for the 
Balmer line and ISO/sub-mm measurements to find values 
several times higher than the UV continuum at all redshifts. 
The rise to = 1 is preserved, arguably the fall off at 2 > 2 
is preserved though given the random errors and the system- 
atics in the dust correction no change for 2 > 1 would also 
be consistent with the data. Whether star-formation peaks 
at intermediate redshift {z = 1-2) or continues to high red- 
shifts (2 > 4) is an important test of heirarchical formation 
scenarios (e.g. Baugh et al. 1998). From the current data the 
question must remain open. 

The agreement between the far-IR/sub-mm measure¬ 
ments and Balmer measurements is particularly impressive 
since both attempt to compensate for dust in different ways. 
Flores et al. hnd an upward dust correction of 2.9 ± 1.3 at 
2 = 1 and extinctions of Av = 0.5-0.9 mags, both consistent 
with our best estimates. It should be noted that indepen¬ 
dent ISO observations of the Hubble Deep Field North by 
Rowan-Robinson et al. (1997) give a conflicting value sev¬ 
eral times higher than Flores et al; however the latter is 
derived from a 19x larger area of sky and so is probably 
better determined. 

An important point is that the dust corrections are not 
many times larger than we have found. Larger corrections 
(e.g. xl5) have been argued for by authors such as Meurer 
et al. (1997) based on amounts of obscuration in powerful 
starburst galaxies locally. However it is not clear that the 
local massive starburst galaxies (e.g. Arp 220) are the right 
analogs of the higher redshift field galaxies. When lower lu¬ 
minosity starburst galaxies are studied locally they prove 
much less obscured (e.g. Buat & Burgarella 1998), these sys¬ 
tems are comparable in extinction to local spirals and the 
sample studied here. We note though if we use the dust cor¬ 
rection from the Calzetti law the Ha point at 2 = 1 moves 
up still further to 9x above the UV point. This is because if 
our galaxies are similar to those studied by Calzetti the HII 
regions suffer additional obscuration relative to the stellar 
continuum. This higher value would then agree with the ISO 
observations of Rowan-Robinson et al.. However the result¬ 
ing high star-formation rates may be problematic in that 
the stellar density produced in the local universe might be 
too high (Madau et al. 1998). 


9 CONCLUSIONS 

From a sample of 13 galaxies observed with CGS4 we have 
performed the hrst Ha measurements of the star-formation 
rate at 2 = 1. We conclude the following: 

(i) The Ha measurements show a star-formation rate at 
least three times as high as that inferred from the 2800A 
continuum luminosity by Madau et al. (1996). 

(ii) The typical dust extinction derived using standard 
extinction laws is only moderate {Av = 0.5-1.0 mags) and 
very similar to that inferred in low redshift held galaxies. 
If there is a large population of obscured star-forming sys¬ 
tems at 2 = 1 they are not common in known samples — 
their rate of occurrence in the CFRS at 2 = 1 must be 
< 10% or they would be detected in our data. This limit is 
similar to the results of Pettini et al. for the 2 > 3 Steidel 


et al. galaxies and consistent with deep sub-mm observations 
which estimate that massive obscured star-forming systems 
make up approximately 10% of high-redshift galaxies (Lilly 
et al. 1998B). 

(iii) A cautionary note is the nature of the dust extinction 
law: if we follow the Calzetti attenuation prescription we 
imply much higher obscuration of the Ha line and a star- 
formation rate at 2 = 1 three times higher still. It is unclear 
however whether such a large correction should be applied 
to the integrated light of all galaxies at high-redshift. This 
issue can only be resolved by further direct measurement of 
the H/?/Ha decrement in these galaxies and the derivation 
of the extinction to the nebular regions independent of the 
stellar UV flux. 

(iv) The mean star-formation rate of a 2 = 1 CFRS 
galaxy is ~ 20-60 Mq yr“^ (for Salpeter IMF and the range 
of dust laws we have studied), which is enough to make such 
a L* galaxy in a few Gyr. This is a factor of several higher 
than ordinary spiral galaxies of comparable luminosity to¬ 
day. 

(v) The large scatter in the distribution of Ha to UV 
light is much better ht by a model in which star-formation 
occurs in episodic bursts with intervals of ~ 0.2-0.3 Gyr and 
of length ~ 0.1-0.2 Gyr. Pure continuous star-formation is 
strongly ruled out, even with variable extinction, as a sole 
explanation. 

(vi) We find qualitative trends for star-forming systems 
to have blue colours and peculiar morphology (especially 
interactions). 

(vii) The star-formation history of the Universe, as in¬ 
ferred from the Balmer lines, is qualitatively similar to that 
inferred from the UV but corrected upwards by factors of at 
least 2-3 at all redshifts. Although the overall form of a rise 
to 2 = 1 is preserved, a compilation of dust-insensitive data 
does not yet demonstrate with certainty whether there is a 
turnover beyond 2 = 1. 

Finally, it is clear that the era of detailed spectral stud¬ 
ies of high-redshift galaxies is upon us, made possible by 
the advent of intermediate to high-resolution spectrographs 
on 4m telescopes. In the next few years with larger samples 
and better near-IR spectrographs on 8 m class telescopes 
(many with the ability to spatially resolve spectra in 2D) it 
will be possible to model the spectral evolution of the high- 
redshift galaxy population in much greater details. In par¬ 
ticular the held will be able to focus much more on detailed 
astrophysics (stellar/dust/gas compositions, star-formation 
rates, dynamics, etc.) rather than simple statistics. 
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Table 1. 

OBSERVED SAMPLE AND FLUX MEASUREMENTS 


CFRS# 

Exposure 

(secs) 

Z 

(mags) 

Iab 

(mags) 

(V - I)ab 

EW [Oil] 
(rest A) 

F(Hq) 

lO'^'^ergs cm”^ 

Jc 

(mags) 

HST Morphology 

00.1579 

1000 

0.811 

22.40 

1.45 

33 ± 6 

0.0 ± 13.9 

19.67 ±0.13 

— 

03.0125 

1000 

0.790 

22.08 

2.08 

19 ± 7 

20.1 ± 7.8 

20.19 ±0.10 

— 

03.0133 

1500 

1.048 

22.45 

1.00 

65 ± 16 

52.9 ±33.9 

> 21.87 

— 

03.0316 

1500 

0.815 

21.98 

2.91 

12 ± 2 

0.0 ±6.8 

> 22.87 

Spiral 

03.0615 

500 

1.048 

22.01 

0.96 

27 ±5 

0.0 ± 24.9 

> 20.87 

— 

03.1534 

2000 

0.798 

22.45 

0.70 

39 ± 8 

68.9 ±23.8 

20.70 ±0.23 

— 

10.1220 

1200 

0.909 

22.36 

0.97 

20 ±4 

75.5 ±16.2 

> 21.69 

Peculiar (merger) 

14.0600 

4000 

1.038 

21.53 

0.69 

27 ±9 

108.2 ± 19.0 

19.08 ±0.19 

Peculiar (close pair) 

14.0818 

1000 

0.899 

21.02 

1.12 

19 ± 2 

0.0 ±30.5 

18.74 ±0.11 

— 

14.1496 

2400 

0.899 

21.80 

1.13 

28 ±6 

0.0 ± 15.7 

19.10 ±0.10 

Compact 

22.0770 

3000 

0.819 

21.78 

1.58 

39 ±5 

42.9 ± 6.4 

20.74 ±0.18 

— 

22.1313 

6000 

0.819 

21.74 

0.84 

40 ±5 

40.5 ± 3.6 

20.36 ±0.08 

Peculiar (merger) 

22.1406 

4000 

0.818 

22.16 

1.16 

55 ±4 

78.3 ± 6.2 

21.24 ±0.20 

Compact 


Table 2. 

DERIVED LUMINOSITIES AND VELOCITY WIDTHS 


CFRS# 

L(Ho) (■) 
(lO^^ergs s“^) 

L{[OII]) (■) 
(lO'^^ergs s ^) 

L(2800A) 

(lO^^ergs s“^ Hz“^) 

FWHM 
rest km/sec 

00.1579 

0.0 ±5.2 

3.8 ±0.5 

18.5 

— 

03.0125 

7.1 ± 2.8 

6.5 ± 1.5 

10.6 

457 

03.0133 

34.7 ± 22.2 

12.9 ±1.8 

52.2 

392 

03.0316 

0.0 ±2.5 

2.4 ±0.4 

3.6 

— 

03.0615 

0.0 ± 16.3 

22.6 ±2.4 

81.3 

— 

03.1534 

24.8 ±8.6 

13.2 ±1.8 

42.7 

283 

10.1220 

36.2 ±7.8 

7.7 ± 1.4 

43.5 

260 

14.0600 

69.5 ± 12.2 

36.4 ±8.3 

154.9 

377 

14.0818 

0.0 ± 14.3 

16.1 ± 1.9 

125.3 

— 

14.1496 

0.0 ±7.4 

14.9 ±1.9 

60.5 

— 

22.0770 

16.3 ±2.4 

13.5 ±1.1 

29.8 

unresolved 

22.1313 

15.4 ± 1.4 

19.5 ±2.5 

70.8 

251 

22.1406 

29.8 ±2.4 

21.8 ±1.2 

33.3 

279 


Notes: 


(i) To correct for dust using the final extinction values derived in Section multiply the above values by the following factors: 
L(Ra): xl.6, L{[OII]): x2.4, L(2800A): x3.1. 

(ii) The instrumental resolution ranges from 70-100 km/s (FWHM) and has been subtracted in quadrature from these values. 
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Table 3. 

CONVERSIONS OF Ha,UV TO STAR-FORMATION RATES 

Luminosities are for 1 Mq yr“^. 


Model 

Z/Zq 

IMF 

L(Ha) 

(10'*^ergs s“^) 

L(1500A) 
(lO^^ergs Hz 

0.1 Gyr 1.0 Gyr 

3.0 Gyr 

L(2800A) 

(lO^^ergs Hz 

0.1 Gyr 1.0 Gyr 

-1) 

3.0 Gyr 

L(2800A)/L(Ha) 
(10-1'^ Hz-1) 

(1 Gyr) 

BC96 {kl96) 

0.02 

SG 

0.68 

4.36 

6.64 

7.19 

2.84 

6.23 

8.63 

9.2 

BC96 {kl96) 

0.20 

SG 

0.58 

3.93 

5.36 

5.40 

2.84 

5.35 

6.58 

9.2 

BC96 (kl96) 

0.40 

SG 

0.50 

3.80 

4.82 

4.83 

2.89 

4.96 

5.80 

9.9 

BC96 (kl96) 

1.00 

SG 

0.40 

3.48 

4.02 

4.02 

2.93 

4.36 

4.81 

10.9 

BC96 (kl96) 

2.50 

SG 

0.28 

3.07 

3.33 

3.33 

2.93 

3.88 

4.10 

13.9 

BC96 (gs95) 

1.00 

SG 

0.61 

3.64 

4.39 

4.40 

2.39 

3.72 

4.13 

6.1 

BC96 (gsHR) 

1.00 

SG 

0.61 

3.64 

4.39 

4.40 

2.39 

3.72 

4.13 

6.1 

PEG (Pad) 

1.00 

SC 

0.41 

4.13 

4.72 

4.72 

3.15 

4.80 

5.39 

11.7 

PEG (Gen) 

1.00 

SC 

0.45 

3.84 

4.42 

4.43 

2.94 

4.57 

5.18 

10.2 

M98 

1.00 

SC 



3.50 



5.10 



BG96 (kl96) 

0.02 

SP 

2.23 

10.21 

12.56 

12.87 

6.29 

9.48 

10.78 

4.2 

BG96 (kl96) 

0.20 

SP 

1.96 

9.27 

10.93 

10.96 

6.32 

8.85 

9.53 

4.5 

BG96 (kl96) 

0.40 

SP 

1.70 

9.22 

10.49 

10.50 

6.59 

8.76 

9.24 

5.1 

BG96 (kl96) 

1.00 

SP 

1.35 

8.61 

9.37 

9.37 

6.85 

8.46 

8.73 

6.3 

BG96 (kl96) 

2.50 

SP 

0.90 

7.62 

7.99 

7.99 

6.87 

7.94 

8.07 

8.8 

BG96 (gs95) 

1.00 

SP 

1.97 

8.73 

9.77 

9.77 

5.50 

7.00 

7.25 

3.6 

BG96 (gsHR) 

1.00 

SP 

1.97 

8.73 

9.77 

9.77 

5.50 

7.00 

7.25 

3.6 

PEG (Pad) 

1.00 

SP 

1.19 

8.64 

9.31 

9.31 

6.16 

7.68 

7.97 

6.5 

PEG (Gen) 

1.00 

SP 

1.28 

8.02 

8.66 

8.66 

5.68 

7.16 

7.46 

5.6 

K83 

1.00 

SP-like 

1.12 








M96 

1.00 

SP 



11.06 



7.04 



M98 

1.00 

SP 

1.41 


8.00 



7.90 




Notes: 


(i) SC and SP are the Scalo (1986) and Salpeter (1955) IMFs 

(ii) K83 is Kennicutt’s (1983) conversion for a Salpeter-like IMF 

(iii) M96 is Madau et al. (1996) values for 0.1-1 Gyr, derived from Bruzual &; Chariot (1993); M98 are the values from Madau et al. 
(1998) 

(iv) BC96 are values derived from Bruzual &: Chariot (1996) models, multi-metallicity ‘kl96’ stellar spectra are based on stellar model 
atmospheres (Lejeune et ai, 1996, Kurucz, 1995) and ‘gs95’ are based on observed Gunn & Stryker (1983) spectra. All use the ‘Padova’ 
stellar evolutionary tracks. 

(v) PEG are values derived from the PEGASE models (Fioc et al. 1997) for the ‘Padova’ and ‘Geneva’ stellar evolutionary tracks 
respectively. 
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Table 4. Best fit model parameters 


Description of model 

Best fit parameters with errors, in the form 
= 2(7 < 1(T < best value < 1(7 < 2(7 

log(L(,est) 

Monte-Carlo deviation 
of L}jest (percentiles) 

Continuous star-formation 

rate 

Ca = 0.61 < 0.65 < 0.70 < - < - 

-38.42 

98.6 

Continuous plus constant 
starbursts of fixed mass 

F = 0.87 < 0.90 < 1.00 < - < - 
Ccr not relevant because = 1.0 

Mf, = 1.23 < 1.31 < 1.39 < 1.48 < 1.57 x 10^ Mq 

T = - < 50.6 < 80.0 < 88.4 < 97.8 Myr 

-36.09 

87.4 

Continuous plus constant 
starbursts of variable mass 

F = 0.93 < 0.96 < 1.00 < - < - 
Ccr not relevant because F^est = 1-0 

= 1.57 < 1.84 < 2.15 < 2.56 < 3.03 x 10^ Mq 
M„ = 0.18 < 0.22 < 0.30 < 0.32 < 0.34 

T = 65.9 < 72.6 < 80.0 < 107.0 < 143.1 Myr 

-35.88 

82.0 

Continuous plus exponential 
starbursts of fixed mass 

F = 0.83 < 0.94 < 1.00 < - < - 
Ca not relevant because F^est = 1-0 

= 4.28 < 4.71 < 5.18 < 5.54 < 5.93 x 10^ Mq 

T = 122.7 < 140.1 < 160.0 < 170.5 < 181.7 Myr 

-36.21 

88.4 

Continuous plus exponential 
starbursts of variable mass 

F = 0.92 < 0.96 < 1.00 < - < - 
Ca not relevant because F^est = 1-0 

= 2.67 < 2.91 < 3.16 < 3.37 < 3.60 x 10^ Mq 
Ma = - < - < 0.10 < 0.13 < 0.17 

T = 65.0 < 72.1 < 80.0 < 106.3 < 141.4 Myr 

-36.04 

88.0 


Notes: 

(i) Parameter ranges ‘Icr’ and ‘2cr’ correspond to Alog>C = —0.5, —1.0. 

(ii) Ma is expressed as a fraction of M^. 
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